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Abstract 

\Q . 

, We provide numerical evidence of global diffusion occurring in slightly perturbed integrable 

Hamiltonian systems and symplectic maps. We show that even if a system is sufficiently close to be 
integrable, global diffusion occurs on a set with peculiar topology, the so-called Arnold web, and is 
, qualitatively different from Chirikov diffusion, occurring in more perturbed systems. 

d ' 

1 Introduction 

The characterization of mechanisms for diffusion of orbits in quasi-integrable Hamiltonian systems and 
^ symplectic maps is a relevant topic for many fields of physics, such as celestial mechanics, dynamical 
i/-) astronomy, statistical physics, plasma physics and particle accelerators. Any dynamical state of an inte- 
grable system can be characterized with a set of action-angle conjugate variables (Ji, . . . , J n , <f>i, . . . , <p n ), 
^ ; where the actions Ij are constants of motion, while the angles 4>j simply change linearly with time. The 
properties of the system which are relevant for the stability are determined by the actions. 

Small perturbations of integrable systems can produce a slow drift of the initial value of the actions, 
5 '• and after certain time these drifts can cumulate in such a way to drive the system into very different 
physical state. By global drift we mean that orbits explore macroscopic regions of the action-space. If 
p. the drift behaves as a diffusion process (precisely, we mean that the numerically computed mean squared 
distance of Ji, . . . , I n with respect to their initial values grows, on average, linearly with time), we will 
refer to it as global diffusion. 

In 1979 Chirikov [1] described a possible model for global drift valid when the perturbation is greater 
than some critical value. Chirikov's model has so far been successfully used to describe diffusion in 
systems from different fields of physics (see, for example, [2]). One of the reasons of the broad detection 
of the Chirikov's diffusion is that its typical times fall within the simulation abilities of modern computers 
as far back as the seventies. 

For smaller perturbations the systems fall within the range of celebrated perturbation theories such as 
KAM [3,4,5] and Nekhoroshev theorems [6], which leave the possibility for global drift only on a subset 
of the possible dynamical states with peculiar topology, the so-called Arnold web, and force diffusion 
times to be at least exponentially long with an inverse power of the norm of the perturbation. The 
theoretical possibility of global drift in quasi-integrable systems has been first shown in 1964 by Arnold 
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[7] for a specific system, and is commonly called Arnold diffusion. Recently, Arnold's result has been 
much generalized [8]. Being interested to applications to specific systems, and in particular to systems 
of interest for physics, we use a numerical approach which, avoiding theoretical difficulties, measures 
directly the quantitative features of eventual long term diffusion. In 2003 we numerically detected a 
very slow local diffusion confined to the Arnold web [9] in a model perturbed system. In this paper 
we provide numerical evidence both on quasi-integrable Hamiltonian systems and symplectic maps of a 
relevant phenomenon of global diffusion of orbits occurring on the Arnold web. More precisely, we show 
that a set of well chosen initial conditions practically explores the whole web and the process behaves as 
a global diffusion. 

In section 2 we will describe the Arnold web of quasi-integrable Hamiltonian systems and symplectic 
maps and in section 3 we will describe the numerical methods used to detect the Arnold web and the 
global diffusion on it. 

2 Arnold web of quasi-integrable systems 

For definiteness, we refer to the Hamiltonian system with Hamilton function: 

1\ J| 1 

-^ = 77 + 77 + 4 + e/ ' f = 72 \ i 72 \ i 1 a \ i q i — ' W 

2 2 cos(0ij + cos(0 2 ) + cos(0 3 ) + 3 + c 

where e is a small parameter and c > 0, so that the equations of motion for Ii,I 2 -,h € K. and 0i, 2 , 03 € 
8 1 are: Jj = — §^ and fa = ^ for any i — 1, 2, 3. In the integrable case (when e = 0) the actions 
are constants of motion while the angles fa(t) = fa(0) + ht, fa(t) = 02(0) + I 2 t, fa(t) = 03(0) + t 
rotate with frequencies uo\ — h, u 2 — h, ^3 — 1- Therefore, each couple of actions I 1: I 2 characterizes 
an invariant torus T 3 . For any small e different from zero, H t is not expected to be integrable. However, 
if e is sufficiently small, the KAM theorem applies ( H e is real analytic and H is isoenergetically non- 
degenerate): for any invariant torus of the original system with Diophantine non-resonant frequencies 1 
there exists an invariant torus in the perturbed system which is a small deformation of the unperturbed 
one. The complement of the set made of these invariant tori is called Arnold web, and in such a set, in 
principle, the motions can exhibit chaotic diffusion. Because Nekhoroshev theorem also applies to the 
system, any eventual chaotic diffusion will occur on very long times that grow at least exponentially with 
an inverse power of e. Precisely, because H is quasi-convex, the following estimates apply: \I(t) — 
J(0)| < ae 1 / 6 , for any |i| < 6exp(e /e) 1 / 6 with a, b suitable positive constants (see [10] and also [16]). 

To describe the topology of the Arnold web, it is convenient to refer to the subset of the phase space 
determined by the section S = {(I 1: I 2 ) € IR 2 ,0i = 0, i = 1,2,3}. From the KAM theorem, it follows 
that any invariant torus cuts the section S in only one point (Ii, I 2 ). Moreover, the Diophantine condition 
on the frequencies implies that these points are outside a neighborhood of lines fciii + k 2 I 2 + k 3 = 
proportional to / y/\k\ T , for any k — (k 1: k 2 , k 3 ) G Z 3 \0. Therefore, the intersection between the Arnold 
web and the section S consists of all these straight lines with a neighborhood which decreases as \k\ 
increases. Any of these straight lines is called resonance, and the integer \k\ is called resonance order. 
The Arnold web is open, dense, and has a small relative measure (proportional to 7). 

The case of quasi-integrable symplectic maps is analogous. In this paper we consider the following 
quasi-integrable symplectic map: 

0i = 0i + ii , 2 = 02 + I 2 

u>2, ^3 are Diophantine if there exist positive constants 7, r such that | J2i > l/\k\ T , with |fc| = J2i l&*l> f° r a ^ 
k = (ki,k 2 ,k 3 ) e Z 3 \0. The Diophantine condition considered in KAM theorem requires r > 2 (for 3 degrees of freedom 
systems) and 7 which suitably rescales with e. 
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r i = / i + e7rr(0i + /i,02 + /2) , I^h + e^ifa + I^fo + h) (2) 

0<pi 0<p2 

where / = 1/ (cos(0i) + cos(0 2 ) + 2 + c), with c > 0. At small e, the KAM and Nekhoroshev theorems 
([11], [12], [13]) apply to this kind of maps. The resonances of this system are defined by the straight 
lines k^i + k 2 I 2 + 2Trk 3 = 0, with ki, k 2) k 3 £ Z\0, and the topology of the Arnold web on the section 
S = {(h, I 2 ) G K 2 , 4>i = 0, % = 1, 2} is that described for Hamiltonian system (1). 



3 Global Arnold diffusion: numerical results 

A precise numerical detection of the Arnold web is possible with the so-called Fast Lyapunov Indicator 
(hereafter, FLI) method [14]. In [15] we showed that the Arnold web of system (1) indeed corresponds 
to the above theoretical description. In this paper, we are able to choose initial conditions which are 
good candidates to diffuse thanks to the accurate detection of the web provided in [15]. In fact, within 
resonances, both chaotic and regular motions are observed. Of course, regular motions do not diffuse, 
and therefore diffusive orbits must be chosen in the small subset of the Arnold web made of chaotic 
motions. But even within chaotic motions, to have the chance to observe the phenomenon one has to 
restrict the choice of initial conditions to the single order resonances, i.e. to the portion of resonant lines 
which are far from the main crossing between resonances. The reason is that stability times are expected 
to be much longer at a distance of order y/e from these crossings (see, for example, [16]). With this 
selection of initial conditions, we have been able in [9] to show that indeed after times that increase 
faster than power law, initial conditions diffuse along the chaotic border of resonant lines. Though the 
computational time used in [9] was already quite long, due to the slowness of Arnold diffusion, it was 
not sufficient to evidence a global diffusion, i.e. the wandering of motions from one resonance to the 
others. In this paper we provide numerical experiments that thanks to a much longer computational time, 
and proper choice of model parameters, display for the first time such global diffusion. 

The crucial parameters to set in the models are e and the value of constant c appearing in the de- 
nominator of the perturbation. In fact, any Fourier harmonic efk — e J /(</>) exp(— i ^ fcj0)d0 of the 
perturbing function / is proportional to e and decreases (asymptotically) exponentially with \ki\. 
The values of the two parameters must be balanced so that e is smaller than the critical value for Chirikov 
diffusion (determined with numerical methods, see [15], [17]) and the value of c is not too large so that 
harmonics with large order ^ \ ki\ produce measurable effects on the finite time scale of our numeri- 
cal computations. We found suitable values for our experiments e = 0.6, c = 2 for the map (2), and 
e = 0.01, c = 1 for the Hamiltonian system (1). The Arnold webs detected with the FLI method corre- 
sponding to these values are shown in figures 1,2. In the figures, for each point of the action plane (7i, I 2 ) 
we plot a color corresponding to the FLI value according to the color scale reported below the figure. 
The intermediate value of the FLI (orange in the color scale) corresponds to KAM tori. Darker regions 
correspond to initial conditions (ii, I 2 ) on the section S which produce resonant regular motions, while 
yellow regions correspond to initial conditions which produce chaotic motions. Each resonance appears 
as a darker orange region with a yellow border, or as a single yellow line (depending on the resonance 
and on the value of the angles chosen to define the section S; see [15] for more details). Therefore, the 
yellow region on the FLI map corresponds to the chaotic subset of the Arnold web, where an eventual 
Arnold diffusion should be confined. 

Following the explained criterion of choosing initial conditions within resonant chaotic motions far 
from the crossings of resonances we have chosen ten initial conditions near (Ii, I 2 ) = (0.316, 0.146) for 
the Hamiltonian case (see figure 2a), and twenty initial conditions near (Ji, I 2 ) = (1.71, 0.81) for the 
symplectic map (see figure la). Then, we computed numerically the map up to 10 11 iterations, and we 
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integrated numerically with a symplectic integrator the Hamiltonian system up to a t ~ 10 time units. 
The results are reported in figures 1 and 2: on the FLI map of the action plane (ii, I 2 ) we plotted as black 
dots all points of the orbits which have returned after some time on the section S. Of course, because 
computed orbits are discrete we represented points on the double section \<pi\ < 0.005, |0 2 | < 0.005 
for both cases, and moreover 3 = for the Hamiltonian system; reducing the tolerance 0.005 reduces 
only the number of points on the section, but does not change their diffusion properties. Figures la,2a 
show only the location of initial conditions (inside the circles), and figures lb,2b show the result after 
the intermediate times 2.10 9 (10 9 for figure 2b). Figures lc,2c,ld,2d show the result after much longer 
times. To properly display such long term evolutions we needed to use a zoomed out map of the action 
plane. 

In both cases, the orbits filled a macroscopic region of the action plane whose structure is clearly that 
of the Arnold web. The orbits have moved along the single resonances, and avoided the center of the 
main resonance crossings, in agreement with the theoretical results which predict longer stability times 
for motions in these regions. The larger resonances (which correspond to the smallest orders \k\) are 
practically all visited, while this is not the case for the thinest ones (which correspond to the highest 
orders \k\). This is in agreement with the theoretical results of [18], which predict that the speed of 
diffusion on each resonance becomes smaller for resonances of high order. Therefore, the possibility of 
visiting all possible resonances is necessarily limited by finite computational time. 

On average, the drift behaves as a diffusion process. In fact, in the case of the map, the average 
evolution of the squared distance of (7i, I 2 ) from the initial datum, reported in fig. 3, increases almost 
linearly with time, so that we can measure a diffusion coefficient (the average slope of the plot) D ~ 
1.7 10~ 10 (we do not report the computation of the diffusion coefficient in the Hamiltonian case, which 
would require a longer computation time to explore a broader region of action space). 

We remark that this diffusion coefficient characterizes the global diffusion process, while diffusion 
coefficients measured in [9] characterized local diffusion along a specific resonance. 

4 Conclusions 

In this paper we detect a macroscopic diffusion of orbits on the Arnold web of quasi-integrable Hamilto- 
nian systems and symplectic maps. The described diffusion phenomenon is very different from Chirikov 
diffusion, where the overlapping of resonances allows diffusion in macroscopic regions of phase space in 
relatively short time scales and without apparent peculiar topological properties of the stochastic region. 
The study of the long-term evolution in quasi-integrable systems, especially the set up of statistical 
methods, must therefore take into account that, at small perturbations, in a subset of phase space of 
peculiar structure there is an important phenomenon of diffusion of orbits. Indeed, it concerns vari- 
ous problems going from the old question of stability of the Solar System to the modern burden of the 
confinement of particles in accelerators. 
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Figure 1: The four panels correspond to the FLI map of the action plane (/]_, I 2 ) for the map (2), with 
initial condition on the section S (see [15]), with different magnifications. The yellow region corresponds 
to the chaotic part of the Arnold web. Moreover, on panel (a) we mark with a circle the location of the 
twenty initial conditions; on panel (b,c,d) we mark with a black dot all points of the twenty orbits 
which have returned after some time on the section S. We consider 2 10 9 iterations for panel (b); 2 10 10 
iterations for panel (c) and 10 11 iterations for panel (d). 
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Figure 2: The four panels correspond to the FLI map of the action plane (/]_, J 2 ) for the Hamiltonian 
system (1), with initial condition on the section S (see [15]), with different magnifications. The yellow 
region corresponds to the chaotic part of the Arnold web. Moreover, on panel (a) we mark with a circle 
the location of the ten initial conditions; on panel (b,c,d) we mark with a black dot all point of the twenty 
orbits which have returned after some time on the section S. We consider 10 9 iterations for panel (b); 
1.2 10 10 iterations for panel (c) and 1.1 10 11 iterations for panel (d). 
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Figure 3: Average evolution of the squared distance of (Ii,h) from the initial datum for the map, 
measured for the points on the section S. The total computation time t = 10 11 iterations has been 
divided in 10 3 intervals. For each initial condition, and for each interval [(n — 1)10 8 , nlO 8 ], we have 
computed the average of the squared distance of (ii, J 2 ) from the initial datum, taking into account only 
points that in the interval [(n — 1)10 8 , nlO 8 ] are on the section S. Then, we averaged over all particles. 
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